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Abstract 

From Hamilton's principle of stationary action, we derive governing equa- 
tions of two-fluid mixtures and extend the model to the dissipative case without 
chemical reactions. For both conservative and dissipative cases, an algebraic 
identity connecting equations of momentum, mass, energy and entropy is ob- 
tained by extending the Gibbs identity in dynamics. The obtained system is 
hyperbohc for small relative velocity of the phases. 



1 Introduction 

The knowledge of governing equations for fluid mixtures is scientifically and indus- 
trially an important challenge. Many authors derived the governing system by using 
axioms of balance of mass, momentum, energy and second law of thermodynamics 
m. The mixtures were considered as a collection of different media co-existing in 
the physical space. For example, the balance law of momentum is given in the form : 

^ [ PaUadD + [ Pa Uan da = [ Ta n da + [ ha dD, Q = 1, 2 (1) 
Ot Jd JdD JdD JD 

where L> is a fixed volume, dD is its boundary, n is the unit normal to dD, pa 
are the densities of components, Uq are the associated velocities, Uan are the normal 
part of the velocities at dD, Ta are the stress tensors and ha are the volume forces 
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associated with internal forces and interaction between components. The principle of 
material frame-indifference requires that and depend on the thermodynamic 
parameters of the mixture and on the relative velocity w = U2 — ui. To include 
the added mass effect into consideration, one should be supposed that and 
depend also on accelerations of phases. The structure of this dependence which is an 
important source of interaction, is not clear: should it be frame-indifferent or simply 
Galilean invariant? The method of balance laws ([Hi does not give a definite answer 
to this question. 

It exists a different approach based on Hamilton's principle pSl-fTl which is used 
for construction of conservative (non-dissipative) mathematical models of continu- 
ous media with complex internal structure. The terms including interaction between 
different components of the mixtures do not require constitutive postulates difficult to 
interpret experimentally. They come from the direct knowledge of a unique potential 
for the mixture. According to [5 1, we call a homogeneous mixture if each component 
of the mixture occupies the whole volume of the physical space, and a heterogeneous 
mixture if each component occupies only a part of the mixture volume. In this paper, 
we consider only homogeneous binary mixtures, but the method can be extended to 
the case of heterogeneous mixtures [6]. The plan for the article is as follows : 

In section 2, we formulate an extended form of Hamilton's principle of stationary 
action allowing us to obtain the governing equations of motion. The Lagrangian is 
the difference between the kinetic energy depending on the reference frame and a 
thermodynamic potential which is a Galilean invariant. The equations of motion in- 
troduce two new vector fields different from the velocity fields by taking into account 
the relative velocity of the components. They play the same role as the velocity field 
does in the case of a single fluid. 

In section 3, the governing equations are extended to the dissipative case without 
chemical reactions. An algebraic identity connecting equations of momentum, equa- 
tions of mass, energy equation and equations of entropy is obtained. This identity 
can be considered as the dynamic form of the Gibbs identity. 

In section 4, we justify the compatibility of the governing equations with the 
second law of thermodynamics. 

In section 5, Pick's law is derived. We prove that Pick's law is not yet a linear 
phenomenological law but a direct consequence of governing equations and Stokes 
drag force hypothesis. 

In Section 6, we check the properties of hyperbolicity of the governing system 
for small relative velocity of phases. 

As a convention, in the following we shall use asterisk "V to denote conjugate 
mappings or covecfori' (vector lines); subscripts a = 1,2 indicate the parameters 
of the Q*'^ component; the symbol / indicates the identity; V means the gradient 
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operator-line; V* means the gradient operator-column; a* b means the scalar prod- 
uct of vectors a, b (the vector line is multiplied by the vector column); a b* means 
the tensor product of vectors a, b (the vector column is multiplied by the vector 
line); A a means the product of the mapping A by a. vector a ; b* ^4 means the 
covector c* defined by the rule c* = {A* b)*; div A denotes the divergence of a 
Unear transformation A which is a covector defined as follows: for any vector a, 

div {A a) = {div A) a. + tr ( A . 



2 Governing equations in conservative case 

In paper [7], we considered a pure mechanical case (without entropy). Now we con- 
sider the general case. We take the Lagrangian of the binary system in the following 
form: 

^ 1 

L = -paul - Pa^a " W {pi, P2, Si, S2,w) (2) 

where in the whole paper the summation is over fluid components (a = 1, 2) and pa 
are the densities of components, are the specific entropies, Uq are the velocities, 
w = U2 — ui is the relative velocity, Q,a are the external force potentials, is a 
potential per unit volume of the mixture. The dependance of W with respect to the 
relative velocity is analog to take into account the added mass effect in heterogeneous 
two-fluid theory as it was done by Geurst [3J. The fact that W depends on two 
entropies is classically adopted in the literature El 111 El 111 • The potential W is 
related with the internal energy U of the mixture through the transformation 

dW 

U = W - —w (3) 
aw 

so that the total energy of the system is Q : 

2 1 

Let us note that to define the internal energy of one- velocity media it is useful to con- 
sider a moving coordinate system in which the elementary volume of the continuum 
is at rest. The total energy of the continuum with respect to this system is called 
the internal energy of the medium. For a two-velocity medium, there is no reference 
frame in which any motion could be disregarded. This is the reason why the stan- 
dard definition of internal energy is dependent on the relative motion of components. 



3 



The formula ^ implies that the internal energy [/ is a Galilean invariant. The de- 
pendence of U (or W) on w is an important property of multicomponent fluid 
mixtures. 

Let X be the Eulerian coordinates, t be the time, be the Lagrangian coordi- 
nates of each component. The mass and the entropy conservation laws in the Eulerian 
coordinates are: 

d Pet d 

+ div{pa Ua) = , -^(paSa) + div (paSaUa) =0. (4) 

In the Lagrangian coordinates, equations (HJl are equivalent to: 

Pa det Fa = PaO O^a) , Sa = SqO O^a), (5) 

where 

" (6) 



is the deformation gradient at Xq,; PaO O^a) and Sao O^a) do not depend on t. 
The relation between the Eulerian and Lagrangian coordinates is given by the local 
diffeomorphism x = (X^, t), where (pa O^a, t) is the solution of the Cauchy 
problem : 

= u((/)Q,,t), 0a(Xa,O) = Xq 

Let X„ = Tpa (x, t) be its inverse mapping {4>a ° i^a = I)- We define the virtual 
motion of the mixture such that [4. JOJ : 

X — $Q, (XQ,,t,£Q,) , Xq, — ^a (x, t,£Q.) , ^a ^ — I 

*Q(XQ,i,0) = (t)a{y^a,t) , *Q(x,i,0) = iJaix,t), 

where Ea belong to a vicinity of zero. The Lagrangian and Eulerian virtual displace- 
ments are defined respectively as : 

d ^a d^a 

= -7{ (x, t, Eq) I e„=0 , = -T, O^a,^, £a) \ ea=0 ■ (7) 

OEa OEa 

The definitions ([6117]) imply the following relation between (5Xq, and 5q,x IHJ: 

(5qX =-Fa SXa (8) 

The variations of Uq (t,x), pa (i,x) and Sq (t,x) are deduced from © - ([71) and 
from the definition of the Lagrangian coordinates X„ : 

dt ' dt 5t " 
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We obtain in Appendix A the values of 6ua (x, t), 6pa (x, t) and Ssa (x, t), where 
(5/(t,x) means the variation of / when t, x are fixed and diVaiSX.a) means the di- 
vergence with respect to the coordinates X.^. We note that in Q we used different but 
equivalent expressions for these variations. Using the definition (2) of the Lagrangian 
L as a function of p^, Uq, Sq, we introduce the following quantities: 



Rn 



dL 

dpa 



1 2 
- U_ 

2 " 



dW 



dpa 



K* = — 



1 dL 



Po 



Pa dUa 

dL 

dSr, 



(-1)° dW 

Pa dw ' 



(9) 



dW 



We note that a best set of independent variables is : pa, ja = Pa^a, PaSa- However, 
in this case, the corresponding derivation should be given in four-dimensional space 
ifTTl . For the sake of simplicity we use the first set of independent variables. The 
last formula defines the thermodynamic temperature 6a of each component which is 
now a dynamical quantity depending on the relative velocity of the components. 
Let u) = D y. [ti , t2\ be the domain in the four-dimensional space (x, t) and Ua 
be its image in the (Xq, t)-space. Here [ti, is a time interval and is a fixed 
domain. We consider Hamilton's principle in the form. 



5aO, = Sa L duj 







under constraints (111) where 5aCL are variations of a associated with the variation of 

da , 



= 'J'a(x, t^£a)- It means that 5aa 



0. We have to emphasis on the 



fact that the domain D is fixed in the physical space. This particularity is related to 
the impossibility to have material volume in general motion of the mixture. Taking 
into account formulae ([8]), variations in Appendix A and definitions (O, we get 



Sad = J (^RaSpa + Pa K* (Ju^ - PaOa^Sa 
J (^Ra diVaiPaoS^a) " PaO K* Fa ^(5Xq,) - PaO 9. 



dio 
9X„ 



(5X„ I du>, 



In the last expression all quantities are considered as functions of (X^, t). Hence, 

n dSaO^ 



id = / 



PaO 



'5Xn 



5Xq, duJo 
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+ / DiVaiPaO G) dlVa = 



where G = {Ra (5Xq, — K* 5Xq) and Diva is the divergence operator in the 
4-dimensional space uJa associated with (Xcj,t). All the functions are assumed to 
be smooth enough in the domain uja and 5Xq = on duja ■ Then, we obtain the 
equations of motion for each component in Lagrangian coordinates: 



Taking into account the identity — — Fq, = 0, we rewrite (ITOl ) in Eulerian 

at ax 

coordinates in the form : 

dcJ^a ^ ^ dR^ + e — (11) 

dt " 9x (?x ° (?x 

If ui = U2 = u , then Kq = u and (fTOl) is equivalent to 

f + V('. + «) = -*V.. | = | + u-v 

where /i is the enthalpy and is an external potential ifTOl . Conservations of the 
total momentum and the total energy are a consequence of the governing equations 
©, (HB: 

X! ^ ( (2"" + ) + + divi^PaVia (K^ - ) - Pa—^ = 

(13) 

The covector K* is an essential quantity; indeed, Pa^a (but not PqUq ) is the im- 
pulse for the a*^ component of the mixture. Morever, for adiabatic motions, the defi- 
nition of potential flows for two-component mixtures is associated with rot Kq, = 
and not with rot = 0. For potential motion, equation (fTTI) yields additional 
conservation laws (T) [TTl : 

(9K* 

^ + V (K* - R^) = 

In the particular case of bubbly liquids, Geurst was the first to carry out a term analog 
to Kq ||3l. Moreover, if rotKo, / 0, the system of governing equations is not 
conservative in terms of Kq , pa and Sa (the number of conservation laws admitted 
by the system is less than the number of unknown variables); but, nevertheless the 
system can be rewritten in conservative form if we add the gradient tensor Fa as 
unknown variable Q. 
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3 Governing equations in the dissipative case and dynamic 
Gibbs identity 

The conservative fluid mixture model presented in section 2 is relevant to the first 
gradient theory fVT\: the forces applied to the continuous medium are divided into 
volume forces and surface forces. In fluid mixture flows it is reasonable to neglect 
the surface friction forces compared to galilean invariant algebraic volume forces. 
The virtual work of dissipative forces applied to the a*^ component is in the 
form 5Ta = f* (5aX. For the same virtual displacement of two components, 

2 

5x = 5qXi = (^qX2 , the total virtual work of dissipative forces is 6T = ^ 5x. 

a=l 
2 

For a solid displacement, the work 5T is equal to zero and consequently ^ = 

a=l 

0. We specify later the behavior of forces f* . Let us introduce the quantities 
Ma, Ba-, S and E such that: 

Ba = + div{pana) 

2 ^ ^ 

S = ^ Pa + "a 

a=l 

We prove in Appendix B the following property: 

Theorem: For any motion of the mixture, we have the identity 

2 



E -Y. A/>a + (K^ - ii,) - 5 = 

This relation is the most general expression of the Gibbs identity in dynamics. Anal- 
ogous identities were obtained earlier for thermocapillary mixtures lH and bubbly 
liquids |13| . For each component of the mixture, equation of momentum and equa- 
tion of mass are in the form 

= 0, B« = (14) 
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The Gibbs identity implies S = E. Hence, the equation of the entropy S = 0, 

E + f>a= (15) 

is equivalent to the equation of the energy E = 0. We note also that the equations 
M* = and Ba = imply conservation of the total momentum of the mixture 
M = (see ^) 



4 The second law of thermodynamics 

In conservative case, the system with two different entropies is closed by dUl. In 
dissipative case we need additional arguments to obtain equations for each entropy 
Sa which could replace equations (011. We take these equations in the form : 

PaOa^ + fa* (u„ - u) + = (16) 

at 

2 

which must be compatible with (fTSl ). Here pu = ^ Pa'^a is the total momentum, 

a=l 

2 2 

p E Pa and E Qa = 0. The last relation means that we have only internal heat 
exchanges between components. Consequently, if 

j: ^(u.-u) + ^ < 0, (17) 

a=l ° " 

we obtain the entropy inequality HI |9l 

Let us note that relation ( fTTl ) is verified if 

fl = k , f2 - fl , A: > 

and 

/I In 
Ql = - TT)' ^1 = ~^2, K > 

The fact that the inverse temperatures (coldness) appear in the closure relations play 
an important role in other applications lfT4ll . 
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5 Fick's law as a consequence of the governing equations 



The governing equations for each component are : 

M„ = + P«K„^ - - p^e^— - = 

For slow isothermal motions (^i = 62 = 9q = const), we can rewrite these equations 
in the following approximate form: 

or 

M* ~ « - f* - n 

where pa = — ^0 ■Sa is the chemical potential for the a*^ phase. Considering 

dpa 

the difference M2 — M* we obtain : 

= !i - !l ^ ^ (18) 

P2 Pi PI/O2 



where f * = — f * and p = p2 — Pi- Equation (1181 ) is the general form of Fick's 
law. So, Fick's law is not a linear phenomenological law but a direct consequence of 
equations of motion and the Stokes drag hypothesis which was previously noticed by 
Bowen with an other model O. 



6 Hyperbolicity of the two-fluid mixture model 

The hyperbolicity of governing equations is very important because it implies well- 
posedness of the Cauchy problem. In mechanical case (when we neglect the equa- 
tions of entropies) we are back to our previous study EITT]]. The only difference is 
the right-hand side algebraic terms f^^ due to the Stokes-like drag forces. Obviously, 
they do not affect the hyperbolicity analysis. The potential W is then a function of 
pi,P2 and w = The Lagrangian is (with ^1^ = 0) 

^ 1 

L=Y1 1; Pc^l-W {pi ,P2 ,w) 



We gave a sufficient condition of the hyperbolicity of system (01, (fT2l) in the multi- 
dimensional irrotational case where rot Kq = 0. Recall the main results we 
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obtained in this case: after a change of variables, system ((l4l).([T2l)) takes the form 



1(f) --(^(E^.J.))^"' 

1(f)- 

where 

G{ai, (T2, jl, j2) L{pi, P2, jl, j2) - X! Pa, with CTq, = - — 

a=l 

The function G is a partial Legendre transformation of L [pi, p2, ii, '^2) with 
respect to the variables pa '■ 

dG_ _ _ dG _ 

o — Pa, r^. — ^ai 

and the system ( fT9l ). (l20l) can be rewritten in a symmetric form |[T5l[T6l[T7]| 

+ = 0, A = A*,B^ ={Br, i = l, 2, 3 (21) 

where 

U = (CJI, C72, Jl, J2) , A = 

and the matrices B^ can be obtained from ([T9l ). (l20l ). If A is positive definite, 
system (|2TI ) is hyperbolic. We proved in EITTl that the conditions 

dw'^ ' dpi ' dpf dp2 \dp1dp2J 

guarantee the hyperbolicity of our system for small relative velocity of phases. Due 
to relation the inequalities (l22l ) mean the convexity of the internal energy U that 
corresponds to a natural condition of stability. Finally, we established that the sta- 
bility implies the hyperbolicity of the governing equations for small relative velocity 
w, provided that rot Kq, = (condition always fulfilled for one-dimensional flows). 
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Appendix A. 

(9X (?X 

The definition of the Lagrangian coordinates implies ° + ° = 0. 

ot ax 

Taking the derivative with respect to Eq, at zero, we obtain the following equation : 

~~di ~~d^ ^ " Ht "'^ 

Equation (|5]l yields : 

Spa (x, t) det Fa (x, t) + PaS (det Fa) = 6Xa (23) 

Using the Euler-Jacobi identity, 6 {det Fa) = det Fa {x,t) tr ^F^^ SF^ 

and 6 Fa (x, t) = - Fa (x, t) 6F~' (x, t) Fa (x, t) , 6F-' (x, t) = , 

ax 

we deduce : 

6 (det Fa )(x,t) = - det F« tr (sF-^ Fa] = -det F„ trf . Or, 



5 det Fa = -det Fa diVaiS^S.^) (24) 
Substituting (l24l) into ( [231 ) we obtain : 

dpa{-X.,t) = Pa dlVa{dX.a) + oX^ = j—^ 

Pao oXa detra 

5Sa (x, t) = (^Xq, 



12 



Appendix B. 

The proof of the Gibbs identity is obtained by summing the following algebraic iden- 
tities a — / : 
For dissipative terms, 
a. 

ff Ui + fa* U2 - fi* Ui - f2* U2 = 

For the external potentials J7q., 
b. 

^ r, , J- f n \ o o - n 

— paiia + dlV {PailaVia) - Pa - Baila - Pa -q^ = 

For the velocity fields , 

^ Q Pa Ua) + [pa ^a (u^ " ^ )^ 

Let us introduce i = — - — . Then the expression, 

aw 

and the three following identities d — f prove the formula, 
d. 

J^dwdpa^,. (dw \ d (dw\ dw (dp^ \ 

E 5^ -^+drv ^— p. u. j -p^ - —J u.-— ^— + d^v (p. u.) j . 
e. 



En d^a . r. „ doiSa „ 

Pa t>a + Pa Uq — fa = U 

a=l 



§w+|;.,:„(,-ir(£u.)p„u„)-(p„|((-ir£)+p.(-iril^^^ 
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